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K^ Abstract 

^P ' We consider a self-gravitating collisionless gas as described by the 

Vlasov-Poisson or Einstein- Vlasov system or a self-gravitating fluid 
ball as described by the Euler-Poisson or Einstein-Euler system. We 
give a simple proof for the flnite extension of spherically symmetric 

\^ ' equilibria, which covers all these models simultaneously. In the Vlasov 

p^ , case the equilibria are characterized by a local growth condition on the 

microscopic equation of state, i.e., on the dependence of the particle 
distribution on the particle energy, at the cut-off energy Eg, and in 
the Euler case by the corresponding growth condition on the equation 

p\ , of state p = P{p) at p = 0. These purely local conditions are slight 

3 ' generalizations to known such conditions. 

1 Introduction 

In astrophysics, matter distributions which interact by gravity arise on many 
different scales. While the Euler-Poisson system can be used as a sim- 
ple model for a single star, a large ensemble of stars such as a galaxy or 
globular cluster where collisions among the stars are sufficiently rare to be 
neglected is typically modeled by the Vlasov-Poisson system; both systems 
are non-relativistic and possess relativistic counterparts. We refer to [U 13] 
for astrophysical background of these systems. In a well known approach to 



constructing corresponding equilibrium solutions, the system under investi- 
gation is — by a suitable ansatz — reduced to a semi-linear elliptic equation 
for the potential or its relativistic analogue. The crucial question then is, 
under which assumptions on the ansatz the resulting steady state has finite 
mass and compact support, since only such states are of possible interest 
from a physics point of view. In the present paper we give a simple proof 
for these finiteness properties which works for all the indicated models si- 
multaneously and covers (and slightly extends) all those cases known from 
the literature where the assumption is purely local at the cut-off energy or 
at p = respectively. 

In order to be more precise we first consider the models where matter is 
described as a collisionless gas; for the necessary details of what we outline 
below we refer to the next section. In the non-relativistic and time inde- 
pendent case the ensemble of particles (stars) is described by its density on 
phase space, / = f{x,v) > 0, x,f S M^, which obeys the Vlasov-Poisson 
system 

vVJ-VU-V,f = 0, (1.1) 

AC/ = 4ttp, lim U{x) = 0, (1.2) 

|a;|— >-oo 



Pix) = / f{x,v)dv. (1.3) 

Here U = U{x) denotes the gravitational potential and p the spatial mass 
density induced by /; we assume that all the particles have the same mass 
which we normalize to unity. Clearly, the particle energy 

E = E{x, v) = hv\^ + U{x) (1.4) 

satisfies the Vlasov equation (jl.ip , and hence any function of the form 

/ = ^{E) (1.5) 

with a suitable, prescribed function does as well. The time independent 
Vlasov-Poisson system thus is reduced to the semi-linear Poisson equation 

AU = 4tt f (l)(l\v\'^ + U] dv, limC/(x) = 0, (1.6) 

J \2 J |a'|-i-oo 

and the question is under what conditions on (J) the latter equation has solu- 
tions and whether the resulting steady states have finite mass and compact 
support. A necessary condition for the latter is that (p{E) = for E > Eq 
where Eq is a suitable cut-off energy, cf. [251 Thm. 2.1]. 



If instead we describe the matter as an ideal, compressible fluid, then all 
that remains of the Euler equations in the static, time independent case is 
the equation 

Vp + pVU = 0, (1.7) 

where the pressure p depends on the mass density p via an equation of state 

p = P{p), (1.8) 



and the gravitational potential obeys the Poisson equation (jl.2|) . If P is 
strictly increasing on [0, oo[ then (|1.7p and (jl.Sp can be used to express p as 
a function of U, and again the system is reduced to a semi-linear Poisson 
equation. 

The results in |6j imply that solutions to these semi-linear Poisson prob- 
lems which lead to finite mass and compact support must be spherically 
symmetric. Hence we do not loose any relevant equilibria if we make this 
assumption from the start. Under this assumption the characteristic flow of 
the Vlasov equation has the additional invariant 

L:=\xx vl"^, (1.9) 



the modulus of angular momentum squared. We generalize the ansatz (jl.Sp 
to 

/ = (l>{E)L' (1.10) 

which allows for a certain anisotropy in the Vlasov case; here / > —1/2. 

Suppose now that we wish to describe the analogous physical systems 
in a relativistic set-up. On the kinetic level we can consider the so-called 
relativistic Vlasov-Poisson system, where the Vlasov equation is changed to 

V^f-VU-V,f = 0, (1.11) 



Vi + H 



while the Poisson equation (jl.2p and (jl.Sp remain unchanged; like all other 
physical constants the speed of light is normalized to unity. The particle 
energy is redefined as 

E = E{x,v) = y/TTW + Uix), (1.12) 

and the same reduction procedure as outlined above applies. For the static 
Euler equations the velocity field is identically zero, and hence there is no 
difference between the Euler-Poisson and the relativistic Euler-Poisson sys- 
tems here. 



The relativistic Vlasov-Poisson system is neither Gahlei nor Lorentz in- 
variant and is only included here to show that our simple proof covers all 
the models of this type. The genuinely relativistic case has to be modeled 
in the context of general relativity. Assuming spherical symmetry we use 
Schwarzschild coordinates and write the metric in the form 

ds'^ = -e^'^dt^ + e^^dr^ + r^ {dO^ + sin^ V) , 

where Schwarzschild time t coincides with the proper time of an observer 
who is at rest at spatial infinity, r > is the area radius, and the polar 
angles 9 G [0, tt] and (^ G [0, 2tt] coordinatize the orbits of symmetry. The 
static Einstein- Vlasov system takes the form 

,x 



-y===-vj-V^TM'fi'--vj = o, (1.13) 

yi + hr ^ 

e~^^ (2rA' - l) + 1 = SvrrV, (1-14) 

e-2^ (2r/i' + 1) - 1 = Svrr^p, (1.15) 

p{r) = p{x) = / a/1 + |up f{x,v) dv, (1-16) 

p{r)=p{x)= f (^' f{x,v)—t=. (1.17) 

Here x = (r sin cose/?, r sin sin 99, r cos 0) so that r = |x| is the Euclidean 
norm of a; G M^, • denotes the Euclidean scalar product on M'^, / is spherically 
symmetric, i.e., f{x,v) = f{Ax,Av) for all A G S0(3), and ' denotes the 
derivative with respect to r. As to the choice of the momentum variable 
v G M'^ which leads to the above form of the system we refer to |181 |23j . As 
boundary conditions we require asymptotic flatness, i.e., 

lim /i(r) = lim A(r) = 0, (1.18) 

and a regular center, i.e., 

A(0)=0. (1.19) 

For the static Einstein- Vlasov system the particle energy takes the form 

E = E{x, v) = e'^(^) VT+W, (1-20) 

and an ansatz of the form p.lOp again satisfies the corresponding Vlasov 
equation. With this ansatz the quantities p and p defined in ()1.16p and 



(J1.17p become functions of /i. Moreover, e~ can, via (J1.14p and (jl.lOp . be 
expressed in terms of p, and the whole system is reduced to a single equation 
for /x which arises from (J1.15p . 

To complete the set of models which we consider we turn to the Einstein- 
Euler system. Since we consider the static case, the Euler equations reduce 
to the single equation 

Vp+{p + p)Vfi = 0, (1.21) 

where the pressure p depends on p via an equation of state like (jl.8p . The 
field equations (|1.14p . (J1.15p together with the boundary conditions (jl.lSp . 
p.l9p remain unchanged. As in the case of the Euler-Poisson system p 
becomes a function of p, and the system is reduced to a single equation for 
p. 

Up to technical requirements we make the following assumptions. In the 
kinetic case, 

cPiE) > c{Eo - Ef 

for E < Eq close to the cut-off energy Eq, where c > and k < 1 + 3/2. For 
the fluid case we require that P{p) is strictly increasing for p > with 

P'ip) < cp'l^ 

for yO > and small, where < n < 3. In passing we remark that if one 
computes the pressure induced by the ansatz (jl.lOp in the isotropic case 
/ = then it can be written as a function of p which satisfies the fluid case 
assumption with n = k + 3/2, and the restrictions on the growth rates fit. 

We compare our result with known results from the literature. The clas- 
sical example in the context of the Vlasov-Poisson system are the polytropic 
models where 

f{x,v) = {Eo-E)iL'; (1.22) 

the subscript -|- denotes the positive part. The resulting semi-linear Poisson 
equation is the Lane-Emden-Fowler equation. Based on the analysis in [27] 
the corresponding steady states were analyzed in [2j. Here k,l > —1 with 
k + l + 3/2 > 0. Compactly supported steady states are obtained for A; < 3/-|- 
7/2, for k = 3/ + 7/2 the mass is still finite but the support is M^, and for k > 
3/ + 7/2 the mass becomes infinite. In |171l24j extensions of these and related 
results were given for the Einstein- Vlasov system. In [25] it was shown both 
for the Vlasov-Poisson and the Einstein- Vlasov systems that a sufficient 
condition for a compact support is that the ansatz is of the form (|1.22p 
asymptotically for E ^ Eq, but this purely local condition is sufficient only 
if /c < / + 3/2. This result was motivated by and relied on a corresponding 



analysis for the Einstein-Euler system [T3]. In [25] a list of examples from 
the astrophysics literature is given which are covered by such a purely local 
condition, and it is also demonstrated by a suitable counterexample that 
such purely local methods fail for k > I -\- 3/2. The analysis in the present 
paper relies on a purely local characterization of the (microscopic) equation 
of state as well and is subject to the same restriction. Results which are 
not subject to this restriction have been investigated by quite sophisticated 
dynamical systems methods in [5l |9l \TU[ \TT\ . Using a global characterization 
of the ansatz function (j) these results cover polytropes for the full range 
of exponents mentioned above under a size restriction on the initial data. 
Our analysis is much closer to the ones in [14^ [25] , our conditions are less 
restrictive in that only an estimate and not an asymptotic behavior is needed 
at the cut-off, but more important from our point of view is that our proof 
is transparent and short — cf. Section 3 — , and it applies to all the different 
systems specified above. 

There is by now a rich literature on the stability of steady states of 
the Vlasov-Poisson system, cf. [T] [H [131 112] and the references there. It is 
interesting that the character of the stability analysis changes at the thresh- 
old k = 3/2 — k is again a growth rate for the ansatz function and / = 
here — in the sense that below this threshold one can use a reduction pro- 
cedure which gives a stability result simultaneously for the Vlasov-Poisson 
and Euler-Poisson systems while such an approach does not work above this 
threshold [201 [2T1 [22] . We refer to [12] for a complementary instability result 
in the Euler-Poisson case. 

The paper proceeds as follows. In the next section we show in more 
detail how in the static case the models we consider can be reduced to a 
single equation for a suitably defined function y, related to either U or /i. 
The arguments there are known, but we need to put them into a common 
framework. The steady state under consideration has compact support if 
and only if the function y, which starts with a positive value at the origin 
and is decreasing, has a zero. It turns out that in all the cases considered, 
y satisfies an inequality of the form 

y'{r) < -^, (1.23) 

where the mass function m is defined in terms of p by 

r 
2, 



m{r) = 4-7T s p{s) ds, 
Jo 

and p{r) = r'^^g{y{r)) is given in terms of y. In Section 3 we prove that 
under a condition on the behavior of g at y = 0, all functions which satisfy 



(jl.23p have a zero. This rests on two simple observations. Firstly, the mass 
function is increasing since the mass-energy density p is non-negative, and 
secondly, the latter function is in the present context always decreasing, 
up to the possible anisotropy factor r^' in the Vlasov case. In the last 
section we translate the general condition from Section 3 into a condition 
on the microscopic equation of state </> or the macroscopic equation of state 
P respectively and obtain the compact support property for all the models 
considered above. 

2 The basic set-up 

In this section we discuss in more detail how the analysis of steady states 
for the systems which were introduced above can be reduced to that of a 
suitable master equation for the potential or a related quantity. 

2.1 Kinetic models 

2.1.1 The Vlasov-Poisson system 

In the ansatz (jl.lOp a cut-off energy Eq has to be specified which is the 
value of the potential at the boundary of the support of the matter. On the 
other hand we have the standard boundary condition in (jl.2p at infinity, and 
due to spherical symmetry is seems natural to parametrize for a fixed ansatz 
function (j) the solutions by prescribing the value C/(0) of the potential at the 
center. Since this is one free parameter respectively one condition too many 
it is natural to slightly modify the ansatz (ll.lOp . We prescribe a function <l> 
and make the ansatz 

f{x,v) = ^{Eo-E)L' (2.1) 

with / > —1/2, and we look for y = Eq — U with a prescribed value at 
the origin, y(0) = y > 0. Once a solution y with a zero is found we define 
Eq := liyxir^ooyir) and U := Eq — y. In this way the cut-off energy Eq 
is eliminated as a free parameter and becomes part of the solution. The 
following technical assumption on <I> is required for the reduction procedure, 
but it does in general not guarantee the compact support of the resulting 
steady states. 

Assumptions on $. <!• : M — t- [0,cx)[ is measurable, <I*(r/) = for t/ < 0, 
and <I> > a. e. on some interval [0, rji] with r]i > 0. Moreover, there exists 
K > —1 such that for every compact set K cM. there exists a constant C > 
such that 



If we substitute the ansatz (|2.ip into the definition (|1.3p of p we find after a 
short computation that for U{r) < Eq, 



p{r) = cir^^ / ^{Eo - E) {E - C/(r))'+^/2 dE 

Ju{r) 
fEQ-U{r) 

= cir^^ / $(r?) (^0 - U{r) - r?)'+^/^ dr] 



'0 

and p{r) = if U{r) > Eq. Here 



2'+3/2^ / ' ^^ds. 



cr. ~ ... 

Jo Vl - s 
Hence in terms oi y := Eq — U we find that 

p{r) = r^'g{y{r)) (2.2) 

where 

^(y) ,= j Q Pq ^{ri) {y - vY'-^'H^ , y > 0, ^2.3) 

Under the above assumptions on ^ it fohows by Lebesgue's dominated con- 
vergence theorem that g E C(M) Pi C^(]0, oo[) with 

g'{y) = {l + 1/2) Q r c|>(r?) {y - r^y-^/^dr], y > 0, 



and g £ C^(M) if k + ^ + 1/2 > 0. Due to spherical symmetry the semi-linear 
Poisson equation (|1.6p can in terms of y be written as 

l(rV)' = -47rr2'5(y). (2.4) 

In terms of the Cartesian variables we want potentials U £ C'^(U.^) i.e., 
y £ C^(M^). Hence we require that y'(0) = 0, integrate (|2.4p once and have 
to solve the equation 

,, , m(r) 

2/ W = -^ (2.5) 

where 

pr 

m{r) = m{r,y) = 4-7T / s '^ g{y{s))ds. (2-6) 

Jo 

For any y> the equation (|2.5p has a unique solution y £ C^([0, oo[) with 
y(0) = y, and we briefiy review the proof. Firstly, a standard contraction 



argument shows that there is a unique, local solution on some short interval 
[0,S[. This solution extends uniquely to a maximal interval [0, rmax[ where 
by monotonicity, < y{r) < y. If Tmax = oo, we are done, if not, then 
necessarily y(rmax) = 0, and again by monotonicity, y uniquely extends to 
the right via 

y [r) = 2 ' ^ > ^max- 

In addition, y'(0) = 0, and the regularity of g implies that y G C^(M^) as 
desired. 

In the next section we specify a condition on g which guarantees that 
the solution y has a zero at some finite radius R. In the last section we 
translate that condition into one on the ansatz function $, and from y and 
$ we then generate a steady state of the Vlasov-Poisson system which in 
space is supported on the ball with radius R centered at the origin. 

2.1.2 The relativistic Vlasov-Poisson system 



We make the same ansatz (|2.ip as for the Vlasov-Poisson system, with a 
function $ which has the same properties as in l2.1.1l but the particle energy 
is now defined by (J1.12p . We again reduce the full system to the equation 
([23]) . but now iory = Eo-U-l; note that by (fLT2]) . E>1 + U{r). The 
mass function m is defined as in (j2.6p . but in the relation ()2.16p we obtain 
a different form for the function g: 

a(y) ■= / ^i lo ^(^) il + y-r]){il + y- r,f - l)'+'/' dr? , y > 0, ,^ ^. 
^'' \ ,y<0, ^"^ 

where 



q := 2tt / . ds. 

Jo Vl - s 

The function g looks more complicated now, but it has the same properties 
which were stated in the Vlasov-Poisson case, and we arrive at the same 
type of set-up as in I2.1.1[ 

2.1.3 The Einstein- Vlasov system 

First we observe that the unique solution to the field equation ()1.14p which 
satisfies the boundary condition ()1.19p is given by 

e-2^M = l-^^, (2.8) 



where the mass function m is defined as above. This relation defines A only 
as long as the right hand side is positive, a restriction which is due to the 
fact that Schwarzschild coordinates cannot cover regions of spacetime which 
contain a trapped surface. If we eliminate A via (j2.8p and observe that the 
particle energy p.20p depends only on fi, the static Einstein- Vlasov system 
is reduced to a single equation for /i, namely to (J1.15p . 

In order to arrive at a master equation for a suitable quantity y which is 
qualitatively of the same form as before we need to adapt the ansatz to the 
fact that the particle energy (ll.20p is no longer the sum of a kinetic and a 
potential part. Hence we make the ansatz that 

f{x,v) = ^(^l-^^L^, (2.9) 

where <^ has the properties stated in 12.1. 11 We define y := In Eq — fi so that 
e^ = -Eo/e^; notice that the particle energy ()1.20p is always positive so we 
require that Eq > 0. If we substitute the above ansatz into the definitions 
(J1.16p and (J1.17p we obtain the relations 

p{r) = r^^g{y{r)), p{r) = r^^h{y{r)), (2.10) 

where 

giy) ■■= 



cie^y J^-"'" $(7?) (1 - r/)2 {e^y{l - r,)^ - 1)'+^/^ dr],y>0, 

,y<0, 

(2.11) 
with q as in 12.1.21 and 

;,(y) ,= I diev J,'-^-'' $(r?) (e^^l " V? " 1)'^'^' dr^,y>0, ^2.12) 

1 ,y<o, 



with 



di := 2tt / sVl - sds. 
Jo 



'0 

The functions g and h have the same regularity properties as the function 
g in 12.1.11 cf. [251 Lemma 2.2], and the static Einstein- Vlasov system is 
reduced to the equation 



1 f m{r) 

1 — 2m{r)/r \ r^ 



y'(0 = -^^TrrTTVTr ( -^ + 4vrrp(r) ) , (2.13) 



where the mass function m is defined in terms of p as before and p and 
p are given in terms of y by (|2.1Up . For any y > there exists a unique 

10 



solution y G C^([0, oo[) of ()2.13p with y(0) = y which due to the issue of the 
positivity of the denominator is less easy to see, cf. [171 [Mj. 

Once a solution to (j2.13p is obtained, we define yoo '■= linir-xxDyC?^), 
Eq := e^°°, and /i = IhEq — y. Together with the ansatz (|2.9p this yields 
a steady state with all the desired properties, provided y has a zero. It 
turns out that in order to show the latter not the full information of ()2.13p 
is needed, but only the following inequality which immediately follows from 
that equation: 

,, , m(r) 

y\r) < --^. (2.14) 

The difference between the Newtonian, special relativistic, or general rela- 
tivistic cases is then reflected only in the definition of the function g, and 
only an estimate for g at y = which holds in all three cases is needed to 
guarantee a zero for y{r). 

2.2 Fluid models 

2.2.1 The Euler-Poisson system 

We use the static Euler equation (jl.7p together with the equation of state 

p.Sp in order to express p in terms of U. 

Assumptions on P. Let P G C-'^([0,oo[) be such that P' > on ]0, oo[, 

and 

"1 P'{s) 



-ds < oo. 



We define 



fP P'(s) 
Q{P) ■■= / ^ds, p>0, 
Jo s 



so that Q £ C([0,oo[) n C'^{]0,oo[) with Q(0) = and Q'{p) = P'{p)/p for 
p > 0. When written in the radial variable r, (|1.7p reads 

P\p)p' + pU' = Q. (2.15) 

If we divide by p, integrate with respect to r and apply a change of variables 
it turns out that the pair (p, U) satisfies ()2.15p — at least where p{r) > — , 
provided 

Q{p{r)) = c-U{r), r > 0, (2.16) 

with some integration constant c which like the cut-off energy Eq above is 
the value of the potential at the boundary of the matter support. Let 

f°° P'{s) 

ymax := / ds = lim Q{p) g]0,oo]. 

Jo s P^°° 

11 



Then Q : [0,oo[— >■ [0,ymax[ is one-to-one and onto, and we define 

9{y):=[^'y\lll^'-^ (2.17) 

Then g G C(] — oo, ymax[)nC"'^(]0, 2/max[); and writing y = c—U we invert the 
relation (|2.16|) to read as in (|2.2p with / = there. Hence the static Euler- 
Poisson system is reduced to the same equation (j2.5p with mass function 
defined by ([IJ]) with / = and with g defined by (121711 instead of by (fOjl . 
Hence we are in exactly the same situation as in the Vlasov-Poisson case in 
that the crucial question is whether y has a zero R. 

2.2.2 The Einstein-Euler system 

Similarly to l2.2.1l we use the static Euler equation (jl.2ip together with the 

equation of state (jl.Sp in order to express p in terms of fi. 

Assumptions on P. Let P G C^([0, oo[), P > be such that P' > on 

]0, c»[, and 

■' P'{s) 



We define 



s + P{s) 

P' 



ds < oo. 



/■^ P' s) 



SO that Q G C([0,oo[) n C^{]0,oo[) with Q(0) = and Q'{p) = P'{p)/{p + 
Pip)) for p > 0. We rewrite (jl.2ip in the radial variable r, 

P'ip)p' + ip + Pip))f,' = 0. (2.18) 

If we divide hy p + P{p), integrate with respect to r and apply a change of 
variables we see that the pair (p, p) satisfies (|2.18p — at least where p{r) > 
— , provided 

Qip{r)) = c-p{r), r>0, (2.19) 

with some integration constant c. Let 

f°° P'(s) 
2/max := / ——^r^ds = lim Q{p) e]0,oo]. 
Jq s + P{s) P^oo 

Then as before Q : [0,oo[— )• [0,ymax[ is one-to-one and onto, and we define 
g by (j2.17p which has the same regularity properties as before. Writing 
y = c—p we can invert the relation (j2.19p to read as in (|2.2p with I = there. 
Hence the static Einstein-Euler system is reduced to the same equation 
(|2.13p with mass function defined by (j2.6p with / = 0. The issue again is 
whether y has a zero, and this will be determined by the inequality ()2.14p . 

12 



3 The compact-support-Lemnia 

The key to the compact support property for ah the models discussed above 
is the foUowing result. 

Lemma 3.1 Let y € C"'^([0, oo[) with y(0) = y G]0, ymax[ satisfy the estimate 

I, N rair) 
y'{r)< Y on [0,oo[, 



where 



m{r) = m{r,y) := Att s ^ g{y{s))ds, 
Jo 



g G C(] — oo,ymax[) is increasing with g{y) = for y < and g{y) > for 
y > 0, and I > —1/2. Let g satisfy the estimate 

g{y) > cy'"^^ for < y < y* 

with parameters c> 0, y* > 0, and < n < 3 + I. Then the function y has 
a unique zero. 

Proof. Since y is decreasing, the limit yoo := limr_s.oo vif) £ [— oo, oo[ exists, 
and we need to show that yoo < 0. This will imply the existence of a zero of 
y which will be unique by the strict monotonicity of this function. Assume 
that T/oo > 0. Then y(r) > yoo on [0, oo[, and by the monotonicity of g, 

s^^^'ds = ^^giyoo)r^'^^. 

If we put this into the estimate for y' and integrate we obtain the contra- 
diction 

y{i~) ^ y— Cr ~^ — 7- — oo as r — 7- oo; 

C denotes a positive constant which may depend on all the parameters, may 
change from line to line, but never depends on r. 

The argument so far is standard and well known, and the crucial task 
is to derive a contradiction from the remaining possibility that yoo = 0. 
Firstly, we observe that m is increasing in r and positive for r > 0. Hence 

m{r) > m{l) =: rrii > for r > 1, (3-1) 

and 

f°° , f°° ds nil 

y{r) = — I y'{s) ds > nii j — = — for r > 1. (3-2) 



s^ 
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Secondly, since g is increasing and y decreasing, 

mir) > ingiyir)) J' s^+^'ds = ^r2'+3^(y(r)). (3.3) 

Hence 

Att 
y'{r) < --^-— r2'+ig(y(r)), r > 0. 

By a simple change of variables this implies that for all r > 0, 



'i,W9(i() Jo 9(s(s)) "2/ + 3 Jo (2i + 3)(2i + 2) 

Now we take r > sufficiently large so that < y(r) < y*; recall that by 
assumption, y{r) — )• as r — )• oo. Then by the growth assumption on 5, 



We estimate the left hand side from below using ()3.2p , multiply the resulting 
estimate by y(r)-^'+-^ and compute the integral where we need to distinguish 
the cases n + / 7^ 1 and n + l = 1. In the former case we find that 

Ci < j^^^Y^ {{y*f-'-'' - y{rf-'-') y{rf^' + c^yirf^', 

in the latter 

Cl<-lnf4^')y(rf+^ + C2y(rf+^ 
c \y{r)J 

which holds for r sufficiently large with positive constants Ci,C2- Since 

2/ + 2 > and / + 3 — n > and since by assumption y^o = 0, the right hand 

side goes to zero as r goes to infinity which is the desired contradiction. □ 

Remark. The two estimates (13. ip and (13. 3p on which the above proof rests 

are quite obvious from a physics point of view. The mass function m{r) 

is increasing in r since the mass-energy density p is non-negative, and this 

yields p.ip . and /O is, up to a possible anisotropy factor r^' in the Vlasov 

case, a decreasing function, which yields (|3.3p . 



4 Application to the various models 

In this section we apply Lemma 13.11 to the various models. We have to 
check what type of ansatz function $ or equation of state function P leads 
to a relation between the mass-energy density p and the function y with a 
functional dependence g which satisfies the growth condition in that lemma. 
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4.1 Kinetic models 

Let $ satisfy the assumptions stated in 12.1.11 In addition let 

^{r]) >cr]^ for rj e]0,r]o[ (4.1) 

for some parameters c > 0, r/o > 0, and — 1 < A; < / + 3/2. 

4.1.1 The Vlasov-Poisson system 

For < y < r/o the function g defined by (j2.3p satisfies the following estimate; 
C > denotes a constant which can change from line to line and depends 
only on the parameters above: 



9{y)>C j v'{y-vy+'/^dv = Cy'+'/' I r^\l - v/yy+'/' d^ 



Hence g satisfies the assumption in Lemma [3.11 with < n = k + 2>/2 < 3 + / 
by the assumption on k. 

4.1.2 The relativistic Vlasov-Poisson system 

In this case the corresponding function g is defined in ()2.7p . We observe 
that (1 + y — r/)^ — l = (l+y — r/+l)(l + y — 77 — l)>y — r/ and find that 

g{y) >C r v\l + y - 7?) ((1 + y - //)' - 1)'+^/' dr, 
Jo 

>C ri]''{y-r,y+^/^dri 
Jo 

= Cy^+l+3/2 f\k^i-sy+'/^ds. 
Jo 

Again, g satisfies the assumption in Lemma |3. 11 

4.1.3 The Einstein- Vlasov system 

In this case the corresponding function g is defined in (12. lip . We estimate 
analogously to l4.1.2[ and in addition we observe that for y sufficiently small, 
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l-e-y > y/2 and e^ > 1/2. Hence 

10 



g{y)>C r^ V'' {{I - Vf - e-'yy^'^' dv 



fi-e-y 
Jo 

= c(l - e-?^)'^+'+3/2 /' 5^(1 - s)'+i/2 ds > C/+'+3/2 
Jo 

as desired. We collect the results for the kinetic models into a theorem. 

Theorem 4.1 Let <l> satisfy the assumptions stated in \2.1.1\ and ( [^.jp . Then 
for any y > the reduced equation \2. 5|) or 112.13]} has a unique solution 
y G C^{[0, oo[) TO^/i y(0) = y which has a unique zero R > 0. By 

f{x,v) = ^ i y{r) - -\v\^ I \x x v\^^ 
or 

\2l 



fix, v) = ^(^l + y{r) - x/l + bP 



^ ^ \x yc v\ 
or 



f{x,v) = ^(i- e~y^''^ y^TTW) k X 



v\^' 



a static, spherically symmetric solution to the Vlasov-Poisson or relativistic 
Vlasov-Poisson or Einstein- Vlasov system is defined. This solution is com- 
pactly supported, and its spatial support is the ball with radius R centered 
at the origin. The parameter y is related to the potential U or the metric 
quantity fi via y= U{R) — U{0) or y= fj.{R) — /u(0) respectively. Moreover, 

P,peC{R^)ncHBRm- 

4.2 Fluid models 

Let P satisfy the assumptions stated in 12.2.11 or 12.2.21 In addition let 

P'(/o)<cpi/"forpG]0,po[ (4.2) 

for some parameters c > 0, po > 0, and < n < 3. It turns out that in 
checking the condition for the corresponding function g we need not distin- 
guish between the non-relativistic and relativistic cases, since in both cases 

Q{p) < / —^ds < Cp^l-^ 
Jo s 
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for < p < po- Since for positive arguments, g is defined as the inverse func- 
tion to Q this immediately yields the estimate for g required in Lemma |3.H 
and we can sum up our results for the fluid case. 

Theorem 4.2 Let P satisfy the assumptions stated in \2.2.1\ or \2.2.iA and 
^4-^ - Then for any y G]0,ymax[ the reduced equation /i2.5\) or /i2.13]) has a 
unique solution y £ C^([0, oo[) with y(0) = y which has a unique zero R > 0. 
By 

p = g{v), p = P{p) 

with g defined in \2.2.1\ or \2.2.2\ respectively, a static, spherically symmetric 
solution to the Euler-Poisson or Einstein- Euler system is defined. This solu- 
tion is supported in the hall of radius R centered at the origin. The parameter 
y is related to the potential U or the metric quantity /i via y= U{R) — U{0) 
ory= fi{R) - p(0) respectively. Moreover, p e C{R^) D C^{Br{0)). 

4.3 Final remarks 

1. In the kinetic case it is straight forward to extend the above analysis 
to an ansatz of the type 

f{x,v) = ^Eo-E){L-Loy+ 

or its general relativistic analogue, where Lq > and the other pa- 
rameters are as before. Such an ansatz leads to steady states which 
have a vacuum region at the center. This situation was investigated 
in [T9] by a perturbation argument, and the structure of the resulting 
steady states was studied in [1] by numerical means. 

2. The arguments from Lemma 13.11 can easily be applied to show that 
a solution to the equation (|2.4p with data y{r) = y > 0, y'{r) = y' 
prescribed at some radius r > 0, has a zero to the right of f, provided 
y' < 0. The important point is that again y'{r) < —m{r)/r'^ for r >r, 
where m{r) = A-k f° s'^'^ p{s) ds. This extension will be useful in |16j . 

3. In the Euler case there is a size restriction on y, if ymax < oo, i.e., 
if P grows only sub-linearly for large values of p. If the pressure 
is weak in this sense it cannot support an arbitrarily large potential 
difference between the center and the surface of the equilibrium matter 
distribution. For equations of state which typically arise in physics P 
grows super linearly for large values of p so that ymax = oo. To see why 
no such restriction appears in the kinetic case we consider for simplicity 
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the Vlasov-Poisson case with / = 0. Then p = h{y) = h{g ^(p)), i.e., 
P = ho g^^. A simple change of variables shows that in this case 

^, , ff{hog-^y{s) , f9''(p) h'{t) , 

Q(p) = / ^ — ^-^-^ds = / —TT-dt -^ oo as /9 -s> oo 

Jo s Jq g{t) 

because in the isotropic Vlasov case h' is a positive multiple of g, cf. 
Lemma 2.2]. Hence in the Vlasov case ymax = oo. 



4. In the kinetic case, p G C^(M^), provided k + / + 1/2 > 0, and the 
same is true in the fluid case under a suitable assumption on P. 

5. In the kinetic case the restriction / > —1/2 can be relaxed by assuming 
more on <l>. For example, Lemma 13.11 applies to all the polytropes 
(jl.22p in the Vlasov-Poisson case with k,l > —1, k + I + 3/2 > and 
k < I + 3/2, since in that case p{r) = cr'^^y{r)^ . 

Acknowledgment. The results reported here originate in the first author's 
doctoral thesis [15j . 
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